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Spontaneous time-reversal symmetry breaking in the boundary Majorana flat bands 

Yi Li, Da Wang, and Congjun Wu 
Department of Physics, University of California, San Diego, CA 92093 

We study the boundary Majorana modes for the single component p-wave weak topological super- 
conductors or superfluids, which form zero energy flat bands protected by time-reversal symmetry. 
However, due to the divergence of density of states, the band flatness is unstable under spatial 
variations of Cooper pairing phases, which spontaneously breaks time-reversal symmetry and lifts 
the degeneracy. Staggered current loops are generated near the boundary forming non-quantized 
vortices. This effect can appear in both condensed matter and ultra-cold atom systems. 

PACS numbers: 74.45+c,74.81.-g,74.81.Fa 



Recently, Majorana fermions in unconventional super- 
conductors and pairing superfluids have become a re- 
search focus in condensed matter physics [lH2fj|]. They 
appear at boundaries and in vortex cores exhibiting non- 
Abelian statistics which potentially can be used for topo- 
logical quantum computation [J, y, l2l| - |25j . Various the- 
oretical schemes have been proposed for the realization 
and detection of Majorana fermions [9|, [fl UM, I26h29J . 
such as the interface between 3D topological insulators 
and conventional superconductors 9], the proximity ef- 
fect of superconductivity in spin-orbit coupled quantum 
wires [14J Il5| . Cooper pairing with ultra-cold fermions 
with synthetic spin-orbit coupling [29], and p-wave Fesh- 
bach resonances of single-component ultra-cold fermions 
[28j |. Majorana fermions are also proposed in fractional 
quantum Hall states with filling v = | [30l - l32l ]. The sig- 
nature of the zero energy Majorana boundary mode in 
the transport spectra has been observed in the spin-orbit 
coupled quantum wires |16l - l20{ . 



Many unconventional superconductors in 2D and 3D 
exhibit midgap zero energy Andreev modes on the edges 
or surfaces with suitable orientations. This is because 
along the directions of incident and reflected wavevec- 
tors of the Bogoliubov quasi-particles, the pairing gap 
functions have opposite signs [33I - 435J . These unconven- 
tional superconductors can be viewed as weak topologi- 
cal states. They are topologically non-trivial (or trivial) 
along the directions perpendicular to boundaries with (or 
without) zero energy modes, respectively. 

Interactions can significantly change topological prop- 
erties of edge and surface in topological insulators and su- 
perconductors. For non-chiral systems, these boundary 
modes are low energy midgap states, thus they are more 
susceptible to interactions than the gapped bulk states. 
In the helical edge Luttinger liquids of 2D topologi- 
cal insulators, the edge magnetic fluctuations are much 
stronger than those in the bulk, which can lead to sponta- 
neous time-reversal (TR) symmetry breaking and destroy 
edge modes [36H38J. For the zero energy boundary Majo- 
rana modes along the coupled chains of the ^-topological 
superconductors, it has been found that interactions can 
even generate gaps without breaking TR symmetry [39j. 
Similar effects are also found for the helical edge Ma- 
jorana modes of TR invariant topological pairing states 



|40t |41| . In both cases, the new topological class under 
interactions is classified by a Zs periodicity. 

In this article, we consider the most natural interac- 
tion effects in the boundary states in the p-wave weak 
topological superconductors or paired superfluids: the 
coupling between Cooper pairing phases and the zero en- 
ergy Majorana modes. The degeneracy of these boundary 
modes are protected by TR symmetry, but are vulnera- 
ble under spontaneous TR symmetry breaking. Spatial 
variations of phases of Cooper pairing order parameters 
induce bonding among these boundary modes, and thus 
lower the energy. Staggered current loops are generated 
near the boundary, which split the zero energy Majorana 
peaks. Due to the divergence of the surface density of 
states, this TR symmetry breaking mechanism is robust. 

We begin with a heuristic example of two parallel quan- 
tum wires along the cc-direction. Each of them is a ID 
topological superconductor of single component fermion 
with ^-pairing symmetry. As explained in Ref. [5j , zero 
energy Majorana modes exist near the two ends of each 
wire. For example, for the end at x = L of the i-th wire 
(i = 1,2), the operator for the Majorana mode can be 
expressed as 



li 



dx{uo(x)e 



h l {x)+v Q (x)e^ + ^4{x)} 1 (1) 



where 9i is the Cooper pairing phase in the i-th chain. 
We denote A the magnitude of the bulk gap, Ef and kf 
are the Fermi energy and Fermi momentum, respectively. 
In the case of A <C Ef, the zero mode wavefunction is 
approximated as 



u (x) = v {x) « e 
_ k f a 



?» sin k 



f x, 



(2) 



where £ x = -j-^r is the coherence length along the in- 
direction. Now let us connect two right ends with a 
weak link H,j = —t± JV dx [^[(x)^^) + h.c], where 
w is the width of the link. We assume w ^C £ x such 
that the uniform distribution of the pairing phase within 
each wire is not affected by the link. By expressing 

tjji(x) — uo(x)e l ~ +l ~"/i + ... where "..." represents for 
Bogoliubov eigenstates outside the gap A, we arrive at 
two different contributions for the Josephson couplings: 
the usual one through second-order perturbation process, 
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FIG. 1: The fractionalized Josephson coupling between two 
parallel wires through the Majorana zero modes depends on 
sin ^2 instead of cos ^- . Please note the different orientations 
of the x and ^-directions between Fig. [2] and this figure. 



and the fractionalized one through the zero Majorana 
mode as 



E, 



-Jo cos AO — iJi7i72 sin 



A0 



(3) 
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tj_/A, and J x 



where AO = d\ - 

tj - Jl*-al dxu (x)u (x + Ay). 

The fractionalized Josephson coupling of the weak link 
between two wires with the "head-to-tail" configuration 
7, 42] gives rise to 17172 cos A£ ; where 71,2 are Majo- 
rana modes at two sides of the link. In contrast, for two 
parallel wires, the corresponding J\ term in Eq. [3]is pro- 
portional to sin — . This difference can be intuitively 
understood as illustrated in Fig. [T] We rotate the 2nd 
wire by 180°, which corresponds to a phase shift of the 
2nd wire by 7r due to the p x symmetry, and thus A9 is 
shifted to Ad + ir. This changes the dependence of frac- 
tionalized Josephson coupling from cosine to sine in Eq. 
EH As a result, the Majorana fermion bonding strength 
vanishes at AO = which reflects the fact that these Ma- 
jorana zero modes are protected by TR symmetry. The 
bonding strength varies linearly with the phase difference 
as A8 — > 0, thus Eq. [3] is minimized at nonzero values 
of A6q satisfying sin =20 — ± mirx(jj-, 1). The energy is 
gained from the bonding of Majorana fermions. For the 
general case that AOo ^ ±71", time-reversal symmetry is 
spontaneously broken. However, in the case of jj- > 1, 
A#o = ±7r, the Majorana fermion bonding strength is 
maximal. Because of the 2-7T periodicity of phase angles, 
this case is also TR invariant. 

We can further extend Eq. [3] to a group of parallel 
wires whose ends at x — L are weakly connected, which 
forms ID or 2D lattices of Majorana modes coupled to 
the superfluid phases. The effective Hamiltonian is 



H„ 



in] 



-Jo^cos( 



— i J\ y sin ■ 
(W) 



-lilv 



(4) 



where (W) refers to the nearest neighbor bonds. The 
Jo term favors a global phase coherence with a constant 
phase value in each wire. However, this results in a flat 
band of zero energy Majorana fermions which is highly 
unstable. In bipartite lattices, a staggered phase con- 
figuration 8i = {~) 1 0q can minimize the ground state 



energy by generating a uniform bonding strength, where 
the value of 8' can be obtained self-consistently. For a 
2D lattice of surface Majorana fermions, if the superfluid 
phase distribution of 0i forms a vortex, the plaquette 
in which the vortex core is located exhibit a Z2 vortex 
for the Majorana fermion. The motion of the superfluid 
vortex introduces the dynamic Zi flux for the Majorana 
fermions. The study of phase transitions of Eq. 2] will be 
deferred to a later publication. 

Having presented the spontaneous TR symmetry 
breaking of Majorana fermions with weak links, now let 
us consider a 2D superconductor of a single component 
fermion with p x pairing. We use the tight-binding model 
of the following mean-field Bogoliubov-de Gennes (B-de 
G) Hamiltonian with the open boundary condition 

Hm f = ~ Z^2 W xC l Cl +^ + h&i+ey + h.c.) - nc\ci\ 

i 
~ V /Zi A M+e, C ^+e, Ci+h.C.} 



^E A * 



i+e x Ai t i+e^ 



(5) 



where Ai^+g^ = (c-i + g x Ci}. In conventional superconduc- 
tors, the open boundary only affects the magnitude of 
the pairing order parameter not on the phase distribu- 
tion. The p^-type superconductors are weakly topologi- 
cal in 2D, and as is well-known, if the phase distribution 
of A, ,i+e x is uniform such that TR symmetry is main- 
tained, a flat band of Majorana surface modes appear on 
the boundary along the y-direction. 

However, this Majorana band flatness is unstable when 
the spatial variations of pairing phases are considered. It 
spontaneously breaks TR symmetry and generates finite 
bonding strengths among the Majorana modes. Conse- 
quently, the macroscopic degeneracy of the surface zero 
modes are removed and the energy is lowered. To confirm 
this intuitive picture, we perform the self-consistent so- 
lutions to the B-de G equation Eq. [5) The mean-field 
Hamiltonian is diagonalized to obtain the Bogoliubov 
eigenstates whose eigen-operators are defined through 
c * = J2 n U n(i)ln+Vn(iht, where (u n , v n ) T are the eigen- 
vectors with positive energy eigenvalues E n > 0. The gap 
equations for Aj^+g^ read 

Ai.i+g,, = 77 ^ tanh — -^ [u n (i + x)v* n (i) 

- <0 + x)u n {i)), (6) 

The current along the bond between i and i+e a (a = x, y) 

is 

ji,e a = -j-^2 lm [ v n(i + e a )Vn(i)f(E n ) 

n 

+ u n (i + e a )u* n (i)[l-f(E n )}}, (7) 

where f(E) is the Fermi distribution function. 
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FIG. 3: The edge LDOS spectra for the system in Fig. [2] 
6(ui ± E n ) in Eq. [8] are approximated by ^r)/[(ui ± -E TO ) 2 + r/ 2 ] 
with the broadening parameter r\ = 0.04. a), b), c) and d) 
are the LDOS spectra at points 1, 2, 3, and 4 marked in 
Fig. a), respectively, a) LDOS without self-consistency by 
fixing Ai t i+g m at a uniform bulk value, b), c) and d) are self- 
consistent results. The corresponding parameter values are 
presented in the caption of Fig. [2] 



FIG. 2: The current and phase distributions from the self- 
consistent solutions to Eq. [5] with the open boundary condi- 
tion. The system size is L x x L y with L x — 20 and L y — 12, 
and parameter values t x — 1, t y — 1.2, V = 3, and /i = 0. 
Only half of the system along the x-direction is depicted, and 
the distributions in the rest half can be obtained by perform- 
ing the reflection operation, a) The spontaneously generated 
edge current distributions; b) the lengths and directions of 
arrows represent the magnitudes and phase angles #i,;+e x for 
the pairing operators Ai ii+ e x . 



ternal magnetic fields. The size of current loops should be 
at the order of coherence lengths which can be estimated 
as £ a ~ 2t a /A with a = x,y. Currents decay exponen- 
tially along the x-direction at the scale of £ x , and change 
directions along the y-direction at the scale of 2£ y . 

Next we calculate the local density of states (LDOS) 
near the boundary as depicted in Fig. [3] The expression 
for LDOS reads 



L(i,u) 



E 



\u n (i)\ 2 S(to - E n ) + \v n (i)\ 2 S(u + E n ). (8) 



The spontaneous TR symmetry breaking effect ap- 
pears near the boundary normal to the x-direction: cir- 
culation current loops as depicted in Fig. [5] a) . As pre- 
sented in Fig. [2] b) , this edge induces non- uniform phase 
distribution of the pairing order parameter Aj^+g.,. • How- 
ever, it is non-singular which does not exhibit the vortex 
configuration and thus cannot give rise to the current 
loops by itself. Loosely speaking there are two different 
contributions to the current expression of Eq. [7] the 
Bogoliubov modes whose energies outside the bulk gap, 
and the in-gap states of Majorana modes. The former 
contribution can be captured by the Ginzburg-Landau 
formalism as spatial phase variations of A^j+g.,. , while 
the latter arises from Majorana modes cannot. As a re- 
sult, the edge currents do form vortex type configurations 
but actually they are not quantized. The circulation pat- 
tern of these induced edge currents are staggered, which 
is natural, since in the background of superfluidity the 
overall vorticity should be neutral in the absence of ex- 



For comparison, we first present the result without self- 
consistency in Fig. [3] a) by fixing order parameters 
Aij + g, x uniform with the bulk value. The central peak at 
zero energy represents the edge Majorana states. How- 
ever, with self-consistency, spatial variations of Ai ! i +1 s x 
couple to the edge Majorana modes. The zero energy 
peaks are split as depicted in Fig. [3]b), c) and d). Sites 
2, 3 (Fig. [3] b) and c), respectively) are right on the 
boundary normal to x-direction, and thus their coher- 
ence peaks are strongly suppressed. Site 4 (Fig. [3] d) 
is relatively inside, and thus the midgap peaks are sup- 
pressed. Site 2 locates in the middle of current loops, and 
thus the splitting of the zero bias peak is the largest. 

Due to the divergence of density of states at zero en- 
ergy, the above TR symmetry breaking mechanism is 
general. It also applies to the 3D case with the p 2 -pairing. 
Then the xy-surface is non-trivial along which Majorana 
zero energy flat bands appear with the assumption of 
TR symmetry. Again due to the same reasoning, this 
degeneracy will be lifted by spontaneous TR symmetry 



breaking. In this case, we expect that the 2D vortex-anti- 
vortex pattern in Fig. [5] will change to that of 3D closed 
vortex rings, which may further form a lattice structure 
parallel to the xy-suriace. Further study of this lattice 
structure will be deferred to a later publication. More- 
over, the above physics also applies in the continuum not 
just in the tight-binding model. Staggered vortices or 
closed vortex rings will appear near the edge or surfaces 
within the size of coherence lengths, which will also split 
the zero energy peaks of Majorana fermions. 

Next we briefly discuss the boundary Majorana flat 
bands for the TR invariant spinful p^-wave Cooper pair- 
ing. For example, let us consider the pairing order pa- 
rameter in the bulk as Aj,j+ ea . = i\A\<T2& ■ d. With- 
out loss of generality, we can assume that the ci-vector 
lies in the xy-plane with the azimuthal angle <f>, i.e., 
d x + id y = \d\e 1 ^ . This corresponds to that spin-up 
fermions are paired with spin-up, and spin-down ones are 
paired with spin-down, with the relation A-f-j- = |A|e I(?i 
and An ~ \A\e~ 1 ^. Then the boundary Majorana modes 
consist two branches 7-f and 7^ forming Kramer pairs. 
The bonding among Majorana modes can occur within 
each branch of Kramer pair by developing currents as 
explained above. If a residue magnetic interaction ex- 
ists, another way to remove the degeneracy is to develop 
magnetism with the order parameter 27^74,. The spin po- 
larization for such an order is in the a;y-plane. In general, 
these two different bonding mechanisms can coexist. 

The above real p-wave Cooper pairing of single com- 
ponent fermion can be realized in ultra-cold dipolar 
fermionic molecular systems. For example, the 40 K- 87 Rb 
systems have been cooled down below quantum degen- 
eracy [43, H|- Moreover, the chemically stable dipo- 
lar fermion molecules of 23 Na- 40 K have also been laser- 
cooled [45] . If the dipole moments are aligned along the 
z-axis by external electric fields, the interaction between 
two dipole moments exhibits the d r 2_ 3z 2-type anisotropy. 
It means that the interaction is attractive if the displace- 
ment vector between two dipoles is along the z-axis, and 
repulsive if it lies in the xy-plane. Naturally, this leads to 
the p z -type Cooper pairing which slightly hybridizes with 
other odd partial wave channels as predicted in previous 



works [46H48J . For the two-component dipolar fermions, 
the triplet p. Cooper pairing has also been theoretically 
predicted [49, [50] ■ The above predicted spontaneous TR 
symmetry breaking effects will appear on the surfaces 
perpendicular to the z-axis. 

In the condensed matter systems, the boundary Majo- 
rana fermion has been detected in the ID quantum wire 
with SO coupling |16l - [2rj |. A magnetic field is also ap- 
plied in parallel to the wire to lift the degeneracy at zero 
momentum such that the system is effectively one com- 
ponent if the Fermi energy lies inside the gap. The prox- 
imity effect from the superconducting electrodes induces 
superconductivity in the ID wire, which leads to the zero 
energy Majorana ending states. If we put an array of 
these ID wires together and their ends are weakly con- 
nected, these Majorana modes will also spontaneously 
build up bonding strength by developing staggered or- 
bital currents patterns, the consequent splitting of the 
zero bias peaks will be exhibited in the tunneling spec- 
troscopy. 

In summary, we have found that the zero energy 
boundary Majorana flat bands are unstable towards 
spontaneous TR symmetry breaking, which results from 
the coupling between Majorana modes and the Cooper 
pairing phases. The divergence of density of states at 
zero energy leads to spontaneous formation of staggered 
current loops, which generates bondings among the Ma- 
jorana modes and lifts the degeneracy. This is a robust 
mechanism which appears even without introducing new 
instabilities. 

We thank Z. Cai for early collaboration on this work, 
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helpful discussions. Y. L., D. W., and C. W. are sup- 
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Note added Part of this work was presented in the 2011 
APS march meeting [51j. Near the completion of this 



manuscript, we learned the nice work 52j, which also 
studies the spontaneous TR symmetry breaking in Ma- 
jorana flat bands. 
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